Shear viscosity of neutron matter from realistic nucleon-nucleon interactions 
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The calculation of transport properties of Fermi liquids, based on the formalism developed by 
Abrikosov and Khalatnikov, requires the knowledge of the probability of collisions between quasi- 
particles in the vicinity of the Fermi surface. We have carried out a numerical study of the shear 
viscosity of pure neutron matter, whose value plays a pivotal role in determining the stability of ro- 
tating neutron stars, in which these processes are described using a state-of-the-art nucleon-nucleon 
potential model. Within our approach medium modifications of the scattering cross section are con- 
sistently taken into account, through an effective interaction obtained from the matrix elements of 
the bare interaction between correlated states. Inclusion of medium effects leads to a large increase 
of the viscosity at densities larger than ~ 0.1 fm^''. 

PACS numbers: 21.65.-|~f,26.60.-f-c,97.60. Jd 
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Viscosity plays a pivotal role in detcrmiming the sta- 
bility of rotating neutron stars. As Chandrasekhar first 
pointed out [H, emission of gravitational radiation (GR) 
following the excitation of non-radial oscillation modes 
may lead to the instability of rotating stars. While this 
effect would make all perfect fluid rotating stars unsta- 
ble, in presence of viscosity dissipative effects damp the 
oscillations, and may prevent the onset of the instabil- 
ity. As a consequence, a quantitative understanding of 
the viscosity of neutron star matter is required to deter- 
mine whether a mode is stable or unstable (for a recent 
review on neutron star oscillations and instabilities see, 
e.g., Ref. and references therein). 

Early estimates of the shear viscosity coefficients of 
neutron star matter were obtained in the 70s by Flowers 
and Itoh, who used the measured scattering phase shifts 
to estimate the neutron-neutron scattering probability 
[^, Ql- Based on the these results. Cutler and Lindblom 
carried out a systematic study of the effect of the viscosity 
on neutron star oscillations, using a variety of different 
models of equation of state (EOS) of neutron star matter 

i- 

The procedure followed by the authors of Ref. Q, 
while allowing for a quantitative analysis of the damping 
of neutron star oscillations, cannot be regarded as fully 
consistent. Ideally, the calculation of transport proper- 
ties of neutron star matter and the determination of its 
EOS should be carried out using the same dynamical 
model. The work discussed in this paper is aimed at 
making a first step towards this goal. 

We have computed the shear viscosity of pure neu- 
tron matter using a realistic nucleon-nucleon (NN) po- 
tential, the Argonne Wig model previously employed 
to obtain the state-of-the-art EOS of Akmal, Pandhari- 
pande and Ravenhall Q • Within our approach, based on 
the formalism of Correlated-Basis-Function (CBF) per- 
turbation theory [^, 0] , medium modifications of the NN 
scattering cross section are also consistently taken into 



account, through an effective interaction derived from 
the same NN potential. 

The theoretical description of transport properties of 
normal Fermi liquids is based on Landau theory [lo| . 
Working within this framework and including the leading 
term in the low-temperature expansion, Abrikosov and 
Khalatnikov ll| obtained the approximate expression of 
the shear viscosity coefficient 
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where p is the density, vp is the Fermi velocity and m* 
and r denote the quasiparticle effective mass and lifetime, 
respectively. The latter can be written in terms of the 
angle-averaged scattering probability lyV) according to 
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Note that the scattering process involves quasiparticles 
on the Fermi surface. As a consequence, for any given 
density p, the scattering probability only depends on the 
angular variables 9 and (j), the magnitude of all quasi- 
particle momenta being equal to the Fermi momentum 
Pf ~ (Stt^p)^/'^. Finally, the quantity appearing in 
Eq.ijT]) is defined as 
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The exact solution of the equation derived in Ref. [ll| . 
obtained by Brooker and Sykes 13, reads 
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the size of the correction with respect to the result of 
Eq.© being 0.750 < (jy/r^^if) < 0.925. 

Eqs.([T])-([S]) show that the key element in the determi- 
nation of the viscosit y is the in-medium NN scattering 
cross section. In Ref. 



To gauge the model dependence of our results, we have 
replaced the full Argonne ujs potential with its simpli- 
fied form, referred to as Wg which only includes the 
six static operators of Eq. ([5]) and the two spin-orbit op- 



IJ], the relation between NN scat- erators L ■ S (g) [1, (r^ • r,)]. This eight operators are the 

minimal set required to describe NN scattering in 5" and 
P states. The corresponding results, represented by the 
dashed line, show that using the v'^ potential leads to a 
few percent change of rfr"^ over the density range corre- 
sponding to 1/4 < (p/po) < 2, po = 0-16 fm~^ being the 
equilibrium density of symmetric nuclear matter. 



tering in vacuum and in nuclear matter has been analyzed 
under the assumption that the nuclear medium mainly 
affects the flux of incoming particles and the phase space 
available to the final state particles, while leaving the 
transition probability unchanged. Within this picture 
W{9, (j)) can be extracted from the NN scattering cross 
section measured in free space, {da/dQ)vEic, according to 
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where m* is the nucleon effective mass and 6 and (f> are 
related to the kinematical variables in the center of mass 
frame through Ecm = Ppi^ — cos6')/(2m), 9cm = 0- 

The above procedure has been followed in Ref. [l^, 
whose authors have used the available tables of vacuum 
cross sections obtained from partial wave analysis fl^ . 
In order to compare with the results of Ref. [31, we 
have first carried out a calculation of the viscosity us- 
ing Eqs.H])-® and the free space neutron-neutron cross 
section obtained from the Argonne vis potential 
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In the above equation 
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where cr^ and are Pauli matrices acting in spin and 
isospin space, respectively, and 
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The operators corresponding to p = 7, . . . , 14 are associ- 
ated with the non static components of the NN interac- 
tion, while those corresponding to p = 15, . . . , 18 account 
for charge symmetry violations. Being fit to the full Ni- 
jmegcn phase shifts data base, as well as to low energy 
scattering parameters and deuteron properties, the Ar- 
gonne Vis potential provides an accurate description of 
the measured cross sections by construction. 

In Fig. [U we show the quantity rjT'^ as a function of 
density. Our results are represented by the solid line, 
while the dot-dash line corresponds to the results ob- 
tained from Eqs.(43) and (46) of Ref. [15| using the same 
effective masses, computed from the effective interaction 
discussed below. The differences between the two curves 
are hkely to be ascribed to the correction factor of Eq.®, 



not taken into acount by the authors of Ref. [15| , and to 
the extrapolation needed to determine the cross sections 
at small angles within their approach. 
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FIG. 1: (Color online) Neutron matter -qT^ as a function 
of density. Solid line: results obtained from Eqs.lIT])-® us- 
ing the Argonne vis potential and m* computed from to the 
effective interaction described in the text. Dot-dash line: re- 
sults obtained from Eqs.(43) and (46) of Ref. [l5| using the 
same m* . Dashed line: same as the solid line, but with the 
Argonne vis replaced by its reduced form v'g. 

To improve upon the approximation of Eq. ([6]) and in- 
clude the effects of medium- modifications of the NN scat- 
tering amplitude, we have replaced the bare NN potential 
with an effective interaction, derived within the CBF ap- 
proach as discussed in Ref. [Tsj . 

The correlated states of neutron matter are obtained 
from the Fermi gas (FG) states through the transforma- 
tion 
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where the operator F, embodying the correlation struc- 
ture induced by the NN interaction, is written in the 
form 
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S being the symmetrization operator. The two-body cor- 
relation functions fij , whose operatorial structure reflects 
the complexity of the NN potential, can be written in the 
form 
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with the 0,"j given by Eq.®. 

The effective interaction Vcs is defined by the relation 
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where H is the full nuclear hamiltonian and T is the ki- 
netic energy operator. Realistic models of H include, in 
addition to the NN potential Vij , a three-nucleon poten- 
tial Vijk needed to account for the measured binding en- 
ergies of the few-nucleon systems, as well as the empirical 
equilibrium properties of symmetric nuclear matter [lot . 
In this work, we follow the somewhat simplified approach 
originally proposed in Ref. (20| , in which the main effect 
of the three-body force is taken into account through a 
density dependent modification of the intermediate range 
part of Vij. Moreover, in view of the weak model depen- 
dence of rjT"^ (see Fig. [1]), the full ?;i8 potential is replaced 
by its reduced form -Ug, and the contribution of the non 
static components is disregarded 

In order to obtain VeS from Ea. ((T3)) the expectation 
value of H in the correlated ground state is evaluated 
at the two-body level of the cluster expansion . The 
resulting effective interaction reads 
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The radial functions /"("ry) of Ea. ([T^ are solutions of 
a set of Euler-Lagrange equations satisfying the bound- 



ary conditions f {rtj > d) ^ I, 
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= 0, for 
5, 6 (see. 



n = 2,3 and 4. and f"{rij 
e.g., Ref. M)- 

The effective interaction of Eq. p4|) has been tested 
computing the energy per particle of symmetric nuclear 
matter and pure neutron matter in first order perturba- 
tion theory using the FG basis. In Fi g, [^ our results are 
compared to those of Refs. 0] and |2l| . The calcula- 
tions of Ref. 01 (solid lines) have been carried out using 
a variational approach based on the FHNC-SOC formal- 
ism, with a hamiltonian including the Argonne vis NN 
potential and the Urbana IX three-body potential [l^. 
The results of Ref. [2l[ (dashed line of the lower panel) 
have been obtained using the v'g and the same three-body 
potential within the framework of the Auxiliary Field Dif- 
fusion Monte Carlo (AFDMC) approach. The results of 
Fig. [5]show that the effective interaction provides a fairly 
reasonable description of the EOS. 

Note that our approach does not involve adjustable 
parameters. The correlation ranges d and dt have been 
taken from Ref. \2^, while the parameters entering the 



definition of the three-nucleon interaction (TNI) have 
been determined by the authors of Ref. [20| through 
a fit of nuclear matter equilibrium properties. 
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FIG. 2: (Color online) Energy per particle of symmetric nu- 
clear matter (upper panel) and pure neutron matter (lower 
panel). The diamonds represent the results obtained using the 
effective interaction discussed in the text in first order pertur- 
bation theory with the FG basis, whereas the solid lines corre- 
spond to the results of Akmal Pandharipande and Ravenhall 
The dashed line of the lower panel represents the results 
of the AFDMC approach or Ref. [p]. 



Knowing the effective interaction, the in medium scat- 
tering probability can be readily obtained from Fermi's 
golden rule. The corresponding cross section at momen- 
tum transfer q reads 



da 

dh 
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Veff being the Fourier transform of the effective poten- 
tial. The effective mass can also be extracted from the 
quasiparticle energies computed in Hartree-Fock approx- 
imation. 

In Fig. [3] the in- medium neutron- neutron cross section 
at Ecm = 100 MeV obtained from the effective poten- 
tial, with p ^ pq and po/2, is compared to the corre- 
sponding free space result. As expected, screening of the 
bare interaction leads to an appreciable suppression of 
the scattering cross section. 

Replacing the cross section in vacuum with the one de- 
fined in Eq. (|15p . the medium modified scattering proba- 
bility can be readily obtained from Eq. ^ . The resulting 
W{d, (f>) can then be used to calculate rjT^ from Eqs.((T|)- 
®- 

The effect of using the medium modified cross section 
is illustrated in Fig. ID Comparison between the solid 
and dashed lines shows that inclusion of medium modifi- 
cations leads to a large increase of the viscosity, ranging 
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of the effect of viscosity on neutron stars oscillations, we 
emphasize that the work described in this paper should 
be seen as a first step towards the development of a gen- 
eral approach, allowing for a consistent calculation of the 
properties of neutron star matter. 

The authors are grateful to V. Ferrari, for drawing 
their attention to the subject of this paper, and to R. 
Schiavilla, for providing a code for the calculation of the 
NN scattering cros section. Useful discussions with I. 
Bombaci are also gratefully acknowledged. 



FIG. 3: (Color online) Differential neutron-neutron scatter- 
ing cross section at Ecm = 100 MeV, as a function of the 
scattering angle in the center of mass frame. Solid line: cross 
section in vacuum, calculated with the v'g potential. Dot- 
dash line: medium modified cross section obtained from the 
effective interaction described in the text at p = 0.08 fm~^. 
Dashed line: same as the dot-dash line, but for p — 0.16 fm~^. 



between ~ 75% at half nuclear matter density to a factor 
of ~ 6 at p = 2po. Such an increase is likely to pro- 
duce appreciable effects on the damping of neutron star 
oscillations. 
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FIG. 4: (Color online) Solid line: density dependence of rjT^ 
computed using the effective interaction described in the text. 
Dashed line: r/T^ obtained from the free space cross section 
corresponding to the v'g potential. 

In conclusion, wc have computed the shear viscosity of 
pure neutron matter using an effective interaction derived 
from a dynamical model that can also be used to obtain 
the EOS. While our results are interesting in their own 
right, as they can be employed in a quantitative analysis 
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